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Abstract. We show how the quantum analog of the Fokker-Planck equation for 
describing Brownian motion can be obtained as the diffusive limit of the quantum 
linear Boltzmann equation. The latter describes the quantum dynamics of a tracer 
particle in a dilute, ideal gas by means of a translation-covariant master equation. 
We discuss the type of approximations required to obtain the generalized form 
of the Caldeira-Leggett master equation, along with their physical justification. 
Microscopic expressions for the diffusion and relaxation coefficients are obtained 
by analyzing the limiting form of the equation in both the Schrodinger and the 
Heisenberg picture. 

1 Introduction 

1.1 Quantum Brownian Motion 

One of the classic problems in open quantum dynamics is the question of quantum Brownian 
motion, asking how a distinguished 'Brownian' quantum particle experiences friction, diffusion, 
and thermalization due to the interaction with an unobserved surrounding liquid or gas. Start- 
ing with the work of Caldeira and Leggett [1] the bulk of studies on this problem treat the 
environment in a phenomenological way, usually by linearly coupling the Brownian particle 
position to a continuous thermal bath of harmonic oscillators, whose spectral density is then 
chosen as to reproduce the desired relaxation and diffusion constants. 

Using the Feynman- Vernon path integral approach [2] these linear models can even be solved 
exactly for some cases [3l|4j[5l[6] . However, as is well known, these non-Markovian dynamical 
solutions also have limitations. Firstly, they usually have to assume that the Brownian particle 
and the environment are initially in a product state, leading to an unphysical initial transient 
dynamics due to the re-adjustment of the energies once the coupling is switched on [7J. Secondly, 
the generic assumption of a linear coupling with the unbounded position operator, leading to 
spatial correlations over any length scale, can be justified at best for a restricted class of initial 
states. It will be valid if the Brownian particle state is close to a classical state, but may lead to 
unphysically large decoherence rates if the Brownian state is characterized by macroscopically 
large coherence lengths [8]. 

In a more concise sense, one should therefore characterize those situations as generic quan- 
tum Brownian motion (i) where one is interested in timescales such that the Markov assumption 
is permissible and (ii) where the motional states considered are close to a classical state in the 
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sense that the coherence scales are not macroscopic. This is the regime of the Caldeira-Leggett 
master equation for free quantmn Brownian motion [1]. It is obtained from the path integral 
approach in a high-temperature limit and, apart from the temperature T, it contains a friction 
constant as a phenomenological parameter, see Sect. [33] below. As in the corresponding case 
of the classical Kramers equation fSect. [X^ . the diffusion constant Dpp is determined by rj and 
T according to Dpp — rjMkBT, with M the mass of the Brownian particle. This is an instance 
of the fluctuation-dissipation theorem. 

An issue of much debate is the fact that the Caldeira-Leggett master equation is not of 
Lindblad-form and therefore does not preserve the positivity of some initial states. As is well 
known [9] , this can be healed by adding a "position-diffusion" term provided the corresponding 
coefficient D^x satisfies D^-x ^ rj^h^ / (16-Dpp), which in turn leads to "momentum localization". 
Various authors proposed the 'correct' value for Dxx based either on formal arguments, on 
improved evaluations of the high-temperature limit, and, above all, on more specific descriptions 
of the environmental interaction process [7l[9lfT0 l [TT 1 [T2 l[T3lfT4] . In view of this, it seems quite 
plausible that a definite answer concerning the value of Dxx can only be given by accounting 
for the environment in a way which is microscopically more realistic. 

1.2 Microscopically realistic environments 

Given the large body of work on quantum Brownian motion it is perhaps surprising that 
attempts to justify the Caldeira-Leggett master equation by a realistic microscopic description 
of the environment received much less attention. The main reason is the notorious difficulty of 
obtaining realistic effective dynamic equations for the Brownian quantum particle in terms of 
the microscopic properties of a given liquid or gas. Provided such a realistic master equation 
is available in Lindblad form, one has then to take an appropriate limit, which will be called 
diffusive limit" below, to end up with a master equation of the form of Caldeira-Leggett plus 
the "position-diffusion" term. Importantly, the coefficients rj, Dpp, and Dxx are then no longer 
phenomenological constants, but they are determined by the microscopic description of the 
particular environment considered. 

The simplest realistic environment in that sense is clearly given by an ideal gas in a thermal 
state. The gas particles then do not interact with each other, but they influence the Brownian 
particle via two-body forces, which should be taken sufficiently short-ranged to permit a scat- 
tering theory description of the interaction processes. The corresponding effective equation of 
motion for the quantum Brownian particle is called quantum linear Boltzmann equatioi^ Such 
a Lindblad master equation was obtained recently by implementing the Markov assumption 
before performing the partial trace over the gas particles in a non-perturbative calculation [161 
[17] . It is supposed to be valid if the gas is sufficiently dilute to justify both the neglect of 
three-body collisions and the Markov approximation. Moreover, unlike the master equations 
for quantum Brownian motion, it is supposed to be valid even for very non-classical motional 
states, such as the superposition states found in an interferometer. 

The purpose of the present article is to use this quantum linear Boltzmann equation as a 
starting point for obtaining the extended Caldeira-Leggett master equation. We will discuss 
what kind of assumptions and approximations are required to end up in that form and how this 
diffusive limit can be justified by using both the operator and the Wigner-Weyl formulation 
of quantum mechanics in the Schrodinger and the Heisenberg picture. In all cases we will find 
that the friction and diffusion coefficients 77 and Dpp are uniquely specified by combining the 
thermodynamic quantities of the gas with the relevant microscopic properties of its constituent 
particles, namely their mass and their suitably averaged cross section. Moreover, we will see that 
the coefficient of the "position-diffusion" term Dxx, which is required for complete positivity, 
is completely determined by r] and Dpp at its smallest possible value, Dxx — rf'f? j (IGIJpp). 

^ This is the quantum analogue of the classical linear Boltzmann equation for a tracer particle [15j . 
It is important not to confuse this equation, which is non-perturbative, with the linearized quantum 
equation for the reduced single particle gas state of a self-interacting quantum gas. 
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The structure of the article is as follows. In Sect. [2] we briefly review the form of the quantum 
linear Boltzmann equation in operator representation and in the momentum basis. Section [3] 
discusses the diffusive limit of the equation in the Schrodinger picture, both in operator form 
and in the Wigner-Weyl phase space representation. In Sect. |4] we formulate the same limit 
in the Heisenberg picture and discuss the equations of motion for the energy and momentum 
operator. We present our conclusions in Sect. [5] 



2 The quantum linear Boltzmann equation 



Before we present the full form of the quantum linear Boltzmann equation let us briefly 
collect some important steps that lead to its derivation. The idea of using a scattering theory 
formulation for obtaining Markovian master equations goes back to the work by Joos and Zeh 
on coUisional decoherence [18'. However, their master equation, which was later formulated 
in a non-perturbative fashion [19ti20j. cannot describe friction, since the Brownian mass M is 
assumed to be infinitely large compared to the gas mass m, so that energy exchange cannot be 
accommodated. An early proposal for finite mass ratios m/M is the master equation of Diosi 
[12j . which is quite close to the present formulation of the quantum linear Boltzmann equation, 
but, as we will see, it differs in some crucial aspects, such as the inferred value of D^x- A 
perturbative form of the present quantum linear Boltzmann equation was obtained in p!6 l [2T | 
I22j . pointing to a connection with the van Hove relation and the dynamic structure factor of 
the gas. While other important contributions dealt with specific aspects [23*, the present non- 
perturbative form is a quite recent result ^17j, based on a monitoring approach for deriving 
Markovian master equations [Mj . 



2.1 Operator form 



In the following, we will denote the density operator for the motional state of the Brownian 
particle by p, and its momentum operator by P. The quantum linear Boltzmann equation then 
reads 



d _ 1 
dt'^ ~ ift 



2M' 



(1) 



where the mapping C which describes the incoherent effects of the gas environment is given by 
fill 



(2) 



Here the integration is over all momentum transfers Q, and for fixed Q also over the perpen- 
dicular plane = {p e R'^ : p • Q = O}. The Lindblad operators have the form 



Lq,p = e 



_ JQ X/h 



L (p, P; Q) 



(3) 



where X is the position operator of the Brownian particle. The first term in ^ thus effects 
a momentum transfer determined by Q. It is important to stress that the form ([2]) of Cp fits 
the general structure of a translation-covariant and completely positive master equation as 
characterized by Holevo [53] (see [^[?7] for a discussion). 

The function L, which is operator- valued in ([3]), contains all the details of the coUisional 
interaction with the gas. It involves the momentum distribution function p, (p) of the gas, 
its number density ngas and the elastic scattering amplitude / {p^^p^, which determines the 
differential cross section 
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as well as the total cross section 



CTtot (Pj) = J dn \f{p^n,p. 
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where n is a unit vector with dn the associated solid angle element. 

In order to specify L let us denote, for any given momentum exchange Q 7^ 0, the parallel and 
the perpendicular contribution of a vector P by P\\q = [P ■ Q) Q/Q^ and by P = P P\\q, 
respectively. With these definitions the function L is defined by 

„ pnZ^Im 1 / / m\ Q m „ 

^ ^) = iik'~' v-^ + + m) Y + M^ii'^ 

X/ (relip^Q.P^Q) - |,rel(p^g,PxQ) + |) . (4) 
The most natural choice for jj, is of course the Maxwell-Boltzmann distribution 




/^/'(P) = ^575^exp ( I (5) 

with — 2m//3 the most probable momentum at temperature T = 1/ (fce/S). Moreover, in @ 
we have denoted the reduced mass by = mM / (M + m) and relative momenta by 

rel [p, P) p - —P. 

m M 

One limiting form of the quantum linear Boltzmann equation is obtained by replacing the 
scattering amplitudes in Q by their Born approximation. This simplifies the equation consid- 
erably, since the Born amplitude depends only on the difference of the momenta, Jb {Pf,Pi) — 
Ib {Pf ~ Pi: 0)1 which removes the operator- valuedness of the scattering amplitudes. In this 
approximation one has 

^« (- ^'^M"' f + °) 

so that the dp^Q-integration in ^ can be carried out. The resulting equation thus reduces to 
the one proposed in [TBlf^ . 

Another limiting form of the master equation is the case of an infinitely massive Brownian 
particle, m/M 0, where it describes no dissipation, but pure spatial decoherence. As one 
expects, the quantum linear Boltzmann equation reduces in this limit to the proper master 
equation for coUisional decoherence pDll28l[?f] . which was recently tested experimentally 



2.2 Momentum representation 

Although the operator form ([2])- ([5]) of £ will be quite useful for the discussion of the diffusive 
limit below, the physics described by the master equation is more easily understood in momen- 
tum representation. Let us denote matrix elements of p in the basis of improper momentum 
eigenkets as {P\p\P') = p (^P,P'). The incoherent part ^ of the quantum linear Boltzmann 
equation then takes the form 

{P\£p\P') = j dQp{P-Q,P' -Q)Ah^{P,P';Q) 

-p (P, P') \Uq [Mi„ {P + Q,P + Q-Q) + Mi„ (P' + Q,P' + Q; Q)] 
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with the complex function Afin (P, P'; Q) defined by 

Mi„ (P, P'; Q)= f dpL {p, P Q;Q) L* {p, P' Q;Q) . (6) 

One can show that for P = P' this function is equal to the rate density, found in the classical 
linear Boltzmann equation, of the particle ending up with momentum P after a momentum 
gain of Q due to a gas colhsion. In other words, we have Afin {P, P, Q) — M^l (-P; Q) with 

Mi (P; Q)^^ f dpoM (Po) (rel (Po QP), rel {p„P~Q)) 
^ g |^ |rel(po-Q,P)|^-|rel(po,P-Q)|^ ^ _ 

We can equally consider the corresponding classical rate density for the particle with momentum 
P to end up at a different momentum, 

AC (P) J dQMi{P + Q;Q) 

= ^ / dpodQp{po)<J {i-e\ipo,P) + Q,Tc\{p„,P)) 
m, J 

^ ^ |^ |rel(po,P)|^-|rel(po,P) + Q|^ j _ 
This can be used to put the master equation into the shorter form 

{P\£p\P') ^ J<^Q {P, P'; Q)p{P~ Q, P' - Q) 

-1 [A4t iP) + Kit {P')] p {P, P') ■ (9) 



For P ^ P' the function A/in is in general complex and cannot be related to a classical rate, 
but still the analogy to the classical case seems quite intuitive. 



3 Diffusive limit in the Schrodinger picture 

We now want to consider the Brownian motion limit of the quantum linear Boltzmann equation, 
using arguments similar to the treatment that turns the classical linear Boltzmann equation into 
the Fokker-Planck equation ^32j . The situation is actually more complicated in the quantum 
case, since one is dealing with operators whose values can be estimated in a meaningful way 
only when suitable matrix elements are considered. 

We will argue that the quantum counterpart of the classical Fokker-Planck equation can 
be obtained by formally expanding the operators in the collision kernel of the quantum linear 
Boltzmann equation up to second order contributions in the canonically conjugate operators X 
and P. As we shall see, the result is not equivalent to a naive application of the correspondence 
principle on the classical result. Such a procedure would simply lead to the original Caldeira- 
Leggett master equation j which does not guarantee to preserve the positivity of the statistical 
operator. The operator expansion holds under conditions analogous to the classical ones [32] . 
which imply in particular that the statistical operator describes a very massive test particle 
not far from thermal equilibrium, that is to say, close to diagonal in momentum representation. 
More specifically, the off-diagonal elements {P\p\P') may differ significantly from zero only for 
AP := |P — P'l < ^ Mj f5. In the position representation the validity of the expansion requires 
that the statistical operator is only coherent over a length of the order of the thermal wavelength 
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of the test particle, so that {X\p\X') is appreciably different from zero only within a range 
given by the thermal de Broglie wave length Ath = \j2T:h? fi / M , i.e., for AX :— \X — X'\ < 
y/2Trh? P / M . Note that since we require m/M <^ 1 it follows that the Brownian thermal wave 
length is much smaller that the thermal wave length of the gas, Ath ^fi^" — \/2nh'^P/m. 



3.1 Operator formulation 



In order to formulate the Brownian motion limit of the quantum linear Boltzmann equation we 
come back to its explicit expression (Hl-Q and we confine ourselves to the case of a constant 



scattering cross-section \f {pf,p. 
' ir E ' 

Cp — n 



linear Boltzmann equation is given explicitly by 

dQ 



mi 47r 



Q 



o'tot/47r. In this case, the operator form of the quantum 
dp (10) 



JQX/h 



m Q 



2 

m Q m , 
— H 1 

2 M 



M 



P||Q PM" P±Q 



m Q 



M 



\\Q 



iQX/h 



Q -P 



where {•,•} denotes the anti-commutator. Specializing to the case of a Maxwell-Boltzmann 
distribution in the gas ^ we can write 



Cp = Ug 



e'^«>^/''exp 



dQ 



m Q 



i i exp I -/3 




lP2 



\Q 



4M2 



f3 



mQ ■ P 



-iQX/h 



Since we are ultimately interested in an expansion up to second order in P, the contributions 
coming from the terms involving squares of the momentum operator in the exponent will simply 
cancel out. We can replace them by unity, leading to the much simpler expression 



Cp 



m gtot 
mi 47r 



e*Q-^/''exp 



f dQ f ^ f m 

pcxp ( -13 



Q 

mQ ■ P 

AMm^ 



m Q 

2 

mQ ■ P 



,-.Q.x/n _ M cxp -/? 



, mQ ■ P 

2Mm, 



For a small mass ratio m/M ^ 1 this yields finally 



Cp = 



O'tot 



TO 47r 



e^'^-^^/^exp 



Q ■ P 

AM 



pexp 



Q 

2 



(11) 



QP 

AM 



-iQ X/h _ 



exp 



-13 



QP 

2M 



which we will use as the starting point for expanding in X and P. The main requirement for the 
expansion is the assumption that the change in momentum of the Brownian particle is small 
compared to the scales involved in its motional state. More specifically, one has to assume that 
for typical values of the momentum transfer Q the relevant matrix elements of the statistical 
operator vanish unless 

^AX « 1 
n 



7 



and 



M 



AP < 1. 



These conditions are both satisfied if the Brownian state is close to thermal and M m since 
the momentum transfer Q is then typically of the order of the momentum of the colliding gas 
particles — y/2m/ (3. This implies in particular that 



Q 
h 



m _ 1 1^1 
2^ ~^^X^-AX 



and 



M 



' m / /? 



^ < 1 



M MM " \ M ^ AP' 
Expanding the terms in the square brackets of (jlip one thus arrives at 

2 



(^)^P.[P.P]] + ^^[X.{P„P}] 



where a term linear in the momentum transfer has been omitted since it vanishes upon the 
integration in pT|) . It follows that one has to evaluate the integrals 



2M m 4tt 

P rigas CTtot 



— / dp^fi {P±Q + — Qt 



Q 



'^6M m 47r 
The coefficient 77 is given by 



JdQQ J ^ dp HP (^Pj_Q- 



(3 Ugas CTtot 

6M m 47r 



dQQe 



-/3QV(8m) 



dp Hp (p_lq) 



_ 16 ImkBT 

— -^f^gas CTtot 1 



27rM2' 



The final result thus takes the form 
3 



1=1 



(12) 



(13) 



which is an extended version of the Caldeira-Leggett equation. In the original master equation 
the third term is absent, D^x = 0, and 77 is a phenomenological parameter. In contrast, the 
diffusive limit of the quantum linear Boltzmann equation yielded a microscopically defined 
friction constant HI]), and also the diffusion coefficients Dpp and D^x are directly related to 77, 



Dpp = rjMkBT 



(14) 



and 



Dxx = rj 



IQMksT 



UlkBT 



D 



pp- 



(15) 



Equation (|14p is the expected expression of the fluctuation-dissipation relation, while the coef- 
ficient of the "position-diffusion term" [P^, [Pi,p]] has just the minimal value required 
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to ensure the preservation of the positivity of the statistical operator with elapsing time [9]. 
Since this new contribution appears only in the quantum case, it cannot be read out from the 
classical Fokker-Planck equation and had to be fixed on the basis of a microscopic derivation 
at the quantum level, as pointed out in [16]. Note in particular, that Diosi's form of the linear 
quantum Boltzmann equation [12] leads to a different expression of the coefhcient Dxx^ which 
is not just a function of 77 and T. 

Concerning the comparison with the classical Brownian motion, the considered case of a 
constant scattering cross-section applies in the classical formulation to perfectly rigid spheres 
[32] . and in that case the relevant classical scattering cross-section is the geometric one, af^^ = 
-kB? , with R radius of the sphere. We have just shown that the same friction coefficient appears 
for a constant scattering cross-section in the quantum case provided R S> A^^*^ , due to the fact 
that the forward scattering contribution then cancels out in the master-equation [33ll34j . To 
see how the particular choice (fT5|) of D^x relates to the classical description it is now helpful to 
consider Eq. pS]) in the Wigner-Weyl phase space formulation. 



3.2 Phase space description 



In view of the appearance of a quantum mechanically required "position-diffusion" term in the 
expression (fT3|) it is of interest to follow the transition from quantum linear Boltzmann equation 
to quantum Fokker-Planck by means of the Wigner function [35] , which despite the fact that 
it is not a proper probability density allows for a classical phase-space picture of the quantum 
dynamics. As usual, we denote the Wigner function associated to a statistical operator p as 



W{X,P) 



K 



(16) 



Equation ([TT|) then reads 



d_ 

di 



W{X,P) 



gtot m 
^ 47r toJ 

exp ( — /3 



dg 



M {P Q) 



exp 
- exp 



2M 

Q-iP + K/2) 



2M 

Q jP K/2) 
2M 



P Q 
P 



K 

~2 
K 



AK 

{2TThf 



P\P-Q- 



\p\p 



{p + ^\p\p- 



iX-K/h 



2 ' 

^^ 

2 ' 



A closed equation for W (X, P) is obtained by inserting the inverse of 



exp —j3 



,Q {P Q) 



exp ( —f3 



CTtot m 



2M 
Q P 

2M 
dQ 



W{X,P-Q) 



cosh 



Q-Vx]W{X,P) 



exp (-Q ■ Vj 



Q 

cosh 



J3h_ 

4M 



X 



exp 



Q P 

2M 



WiX,P) 
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where cosh denotes the hyperboUc cosine. The equation can now be written more compactly as 

(17) 



d_ 

dt 



exp (-g • Vp) - cosh ( —Q ■ Vx 



where we have used the unitary differential operators exp (— Q • Vp) and exp(/3?i/ (4M) Q ■ Vx) 
effecting a shift in momentum and position, respectively, of the arguments of the Wigner func- 
tion. As one can check, Eq. (jl7p differs from the corresponding classical expression of the linear 
Boltzmann equation [36^ just by the hyperbolic cosine term, which in the classical case is re- 
placed by unity, as one would obtain in the naive classical limit h Q. An expansion up to 
second order of the exponential operators appearing in the kernel of (|17p transforms the integro- 
differential equation for the Wigner function into a partial differential equation. It is the phase 
space representation of Eq. (fT3|) and reads, 

d_ 

dt 



— W (X, P) = yyVp • [PW (X, P)) + DppApW (X, P) + D^^AxW {X, P) 



This is the quantum counterpart of the classical Fokker-Planck equation, for the time evolution 
of the probability density /d (X, P) in phase-space, 

^/ci (X, P) = r^Wp ■ (P/ci (X, P)) + DppApU {X, P) . 

The quantum and the classical phase space equations for Brownian motion thus differ again by 
the appearance of a "position-diffusion" term, which is symmetric with respect to the regular 
diffusion term. 



4 Diffusive limit in the Heisenberg picture 

In the previous paragraphs we understood the quantum linear Boltzmann equation as a mapping 
C acting on the statistical operator, thus working in the Schrodinger picture. In the same spirit 
we considered its diffusive limit, which led to the quantum counterpart of the classical Fokker- 
Planck equation. In the following we will consider its adjoint mapping C* for the time evolution 
of observables. It is defined through the relation 

Tr {kCp) = Tr (p£*A) , 

where the trace operation expresses the duality relation between the space of trace class op- 
erators, which contains the states given by statistical operators p, and its dual, the space of 
bounded operators A characterizing observables. 

Using the form (O-® of C this leads to the explicit identification 



£* A ^ dQ dp [it (p^ p; Q) Q-^Q■Vt^J^^^Q■X/h^ p. 



(18) 



--{it (p,P;Q)L(p,P;Q),A} 



The full differential equation for the time evolution of an Heisenberg operator At is thus given 

by 



dt in 



2M 



dQ I dp 



(p, P; Q) e-'Q-^/'^Ate^^-^/'^L (p, P; Q) 



1 



{Lt(p,P;Q)L(p,P;Q),At} 
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We will take the Heisenberg operator to coincide with the corresponding Schrodinger observable 
at t = 0, i.e., Ao = A. 

An important class of Schrodinger picture observables are those which are only functions 
of the momentum operator, A = A (P). The equation of motion in the Heisenberg picture then 
simplifies considerably. 



dQ dp\L{p,P;Q)\' 



Note that C, and equivalently C* , is covariant under translations [22] in the sense that 



C 



^-zbP/h^^ibP/h 



= e-*-P/'"'/: [p]c*P/^ 



It follows that the algebra generated by the momentum operator is left invariant. Consequently, 
Ao — Aq (P) implies At = At{P) for t > and therefore [At, P] = 0. In particular, recalling the 
definition of the quantum and the classical in-rates in ^ and ([7]), respectively, we find that 
observables given by a function of momentum obey 



dt 



At{P) 



dQMfl (P + Q; Q) e'^^-^/^At (P) e'^'^'^ - At (P) 
j dQMfl (P + Q- Q) [At (P + Q) - At (P)] , 



in strict analogy with the classical formulation. 

We now focus on the time evolution of the expectation values of momentum and kinetic 
energy, A (P) = P and A (P) — P^/ (2M), respectively. We shall first obtain the explicit dynamic 
equations and then the limiting form corresponding to the diffusive limit. Starting from the 
expression for the time evolution of expectation values, 

A(A),, A Tr [Apt) = Tr [ACpt) - Tr {ptC* A) , (19) 

we can now exploit the fact that functions of the momentum operator are mapped by C* to 
functions of the momentum operator. Specifically, it is convenient to first evaluate 

C* [A (P)] - / dQ (P + Q; Q) [A{P + Q)-A (P)] , (20) 

which can be easily dealt with as an equation for (D-numbers, by working in the momentum 
basis. Moreover, assuming that the scattering cross-section is invariant under parity transfor- 
mations the scattering rate ([7]) satisfies 

Mfl (P + Q;P) = Af^l {P^P + Q) (21) 
= M^„' i-P ^ P-Q) 
= Mfl i-P - Q; -Q) , 

so that the operator £* [A(P)] has the same parity as ^(P). According to our convention 
momentum or energy transfers are positive when the test particle gains momentum or energy. 
The fact that Mj^' is positive thus ensures the obvious physical requirement that the change 
of momentum is positive when the momentum transfer is positive, and likewise that a positive 
energy transfer increases the energy. In fact, for A{P) = P the quantity [A {P + Q) — A (P)] is 
simply the momentum transfer in the single collision, while for A{P) = P^/ (2M) the quantity 
[A{P + Q) - A (P)] is the energy transfer E {Q, P) = {P + Qf / (2Af ) - PV (2M) in a single 
collision with momentum gain Q. 
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Let us first write the explicit expression for M^^ {P + Q; P) for a gas described by the 
Maxwell-Boltzmann distribution ([5]): 



M^^ (P + Q; Q) = / dpMM P±Q 



m Q m 



m» 2 



M 



\\Q 



xcr (rcl(p^Q,P±Q) - ^,rel(p^Q,P±Q^ ' ^ 



2 --^z 2 

Inserting this into Eq. (j20p invites the more compact notation 

o-(-P±Q,Q) = J ^dpfip (p^q) <t ^rcl (p^q,Pxq) - y,rcl (PiQ,P±Q 



Q 

2 



(22) 



(23) 



so that we can write 
C* [AiP)] = - 




Pm f dQ ^ 



Q 



a(Pi,Q)exp -P 



8m? 



2M2 



2Mm, 



x[A(P + Q)- A(P) 



(24) 



We now proceed to evaluate this expression explicitly for the case of the momentum and the 
kinetic energy observable, by considering the special case of a constant scattering cross-section, 

~ ftot/47r, as in Sect. [3l In the momentum case, A{P) — P, it is advantageous 

to rescale by the reference value Mvp, where — pp/m — y^2/ {(3m) is the most probable 
velocity of the gas particles. 



Mvp 



Ctot 



/nm'^v^M An 



Q / 

dQ^exp I -/3 



mg2 
8m2 



2M 



Q _ g -mQ ■ P 
2Mm, 



(25) 



Similarly, for the kinetic energy A (P) — P^/ (2Af ) the expression ([24]) leads to 



C* 



2M \Mv[3 J 



''gas 



ctot f dQ 



— exp -/? 



8m1 



mP; 



2M2 



2Mm, 



2M 



PQ 



M 



h. 



(26) 



The evaluation of /i and I2 is most easily done by focussing on their functional expression, thus 
determining them as functions of the dimensionless variable 



U 



V 



Mvp vp 

and carrying out the integration over the scaled momentum transfer K 



K = 



Q 



^From the right hand side of ((25)) we have 



/r(C/)=n,..^^^/di^^cxpl-if^(^ 



47r 0r M 



m^vp 



K 



1 UK 



if2 



(27) 
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and writing 



K = K, 



u 



K 



u 



one finds that only K\\u, i.e., the component parallel to J7, contributes to (|7f|) . We are thus 
left with 



Ii (U) = ngas 



gtot VfJ 
An 0r M 



dK- 



K 



K 



exp -K^ 



1 UK 

2 ^ if2 



Taking now the direction of U as the polar axis and denoting by ^ the cosine of the angle 
between K and U we have to evaluate 



h (U) = n 



O-tot W/3 TO, 



27r 



+ 00 



47r Jo 
CTtot / »7r U 



-n. 



■gas 



47r \ m(3 M U Jq 



dKK^ 

+ 00 



+1 



-{K/2+U(,)^ 



(28) 



-1 



Ctot 



erf ( y + [/ 



erf ( y - (7 



C/2 



■ sinh([/X) 



47r V m(3 M 



[1 + 2C/2] e-"' - [1 - 4C/2 - AU'] 



Here, erf(a;) = 27r~5 JJ^ exp (— dt denotes the error function. The expression (pS)) can be 
further simplified and expressed more compactly by means of the confluent hypergeometric 
function iFi, in particular by using the known expression of iFi (a, 7; z) for indexes a — —1/2 
and 7 = 5/2 [37], 



This way ([28)1 finally becomes 



gtot 16 / Stt to,, 



1 5 

2' 2' 



(29) 



(30) 



with iFi (—5, |; —U^) a positive, monotonically increasing function. 

For the rescaled kinetic energy one has to consider the r.h.s. of (p6)) . so that the function I2 
is given by 



2n„ 



O"tot V0 TO, 

' 47r 



difexp -if 2 



1 UK 

2 ^~?f2- 



1 TO, 

2 M 



K 



U K 

K 



As before, we take the direction of U as polar axis and denote by ^ the cosine of the angle 
between K and U . This way, 

+ 00 /. + ! 



2n. 
-2n 



Q-tot / Stt TO, 
^''^ 47r V TO/3 M Jo 
Ctot / Stt to 



dKK^ 



d^e 



-{K/2+uiy 



1 TO, 

2 M 



1 



47r \j m(3 M Jq 



TO, 



M J AU 



-1 

difif2e-^'/4 

erf ( ^ + C/ I - erf ( ^ - [/ 



if 



-2n„ 



CTtot / Stt to, 
' Att \l mf3 M 



[1 + 2C/2] e"^' - [1 - - 4J72] 



/¥crf(C/) 
2 [T" 



sinh(t/is:) 



TO 

M 



f(^[5 + 2{/2] e-^' + [3 + 12C/2 + 4C/4] ^ | 
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Also this expression can be written much more compactly by exploiting again (|29p . as well as 
the known expression of iFi (—5, |; —U^), as given in [37], and the functional relations between 
confluent hypergeometric functions with different indexes [38]. This yields 



1 

and therefore 



j2\ _ o f^tot 16 /Stt r 2 J, / 1 5 2\ 3 m, ^ 3 3 



(32) 

with both iFi (—5, |; — ?7^) and iFi (— |, |; — C/^) positive and monotonically increasing func- 
tions. 

Finally, coming back to the operator expression with the original variables l|25)) - ([26)) the 
time evolution of the momentum expectation values satisfies, according to p9|) . 




Similarly, we have for the kinetic energy E = P / (2 A/) 



d , , _ _ CTtot 16 / Stt m 



E>P.=-2n,.3^^^/^^ (34) 



/rni 15 ^^m \ ^ 3 ^ f 33 „,-m 
" I -2' 2 = -^^M ) ^ - W^^^^ (-2' r ~^^M 



ipf 



^From this it is immediately clear that in general there is no closed evolution equation for 
either the first or the second moment of the momentum operator, (P)pj or (P^)pt, since due 
to the presence of the confluent hypergeometric functions moments of arbitrary high order are 
involved in the equation. 

However, the equations do get closed in the limit of a very massive test particle close to 
thermal equilibrium, since in this case the velocity V ^ P /M of the test particle is much smaller 
than the typical velocity vp of the gas particles. Thus using V /vp ^ 1 and m/M <C 1 in ([33]) 
and (|34p the confluent hypergeometric functions are replaced by unity, since li^i (a,7;0) = 1, 
and the reduced mass to, is replaced by m. This leads to 



_d CTtot _32_ r^J^/p^ 

dt^ - Att 3mV P ^ 

which describes velocity-proportional friction with coefficient 77 leading to the expected expo- 
nential relaxation to a mean momentum equal to zero. Similarly, 



ii/F\ - _9 CTtot _32_ I 2nm f _ _3_ 
dV^"'^ ^s"^4^3AfV P V^"' 2/3 



-2r, ( (E>p, - — 



3 

2f3 



shows that the mean kinetic energy relaxes exponentially to the equipartition value i/ik^T. 
The relaxation rate 



16 mk^T 
^^yWtot^^, 

is equal to the result (fT2]) obtained in Sect. [3] for the considered case of a constant scattering 
cross-section. 
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5 Conclusions 

In conclusion, we discussed how the "diffusive Umit" of the quantum version of the hnear 
Boltzmann equation yields the master equation for quantum Brownian motion, and provides a 
microscopic formulation of the relaxation and diffusion constants. In particular, we saw that this 
procedure leads naturally to the minimal extension required to turn the Caldeira-Leggett master 
equation into Lindblad form. The approximations invoked in this limit could be physically 
justified by using both the operator and the Wigner-Weyl formulation of quantum mechanics 
in the Schrodinger picture, while the Heisenberg picture provided a dynamic description of the 
relaxation behaviour. Still, a mathematically more rigorous treatment of the diffusive limit is 
clearly desirable, together with an extension of the obtained results to the case of an arbitrary 
scattering cross-section. 
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